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SUMMARY 

During  the  second  year  of  this  effort,  we  focused  on  introducing  the  effect  of  combustion 
in  the  computation  of  the  spatially  developing  shear  layer,  and  continu^  our  analysis  of  the  effect 
of  density  variation  and  upstream  forcing  on  the  growth  of  the  mixing  zone  of  the  shear  layer.  We 
have  developed  a  flame  sheet  model  for  the  simulation  of  combustion  at  high  Damkohler  numbers 
where  the  application  of  the  transport  element  method  proves  to  be  rather  expensive.  The  model 
uses  the  instantaneous  local  strain  rate  as  an  input  from  the  flow  computations  and,  by  integrating  a 
one-dimensional  equation,  computes  the  rate  of  burning  within  each  flamelet  in  the  domain.  The 
coupling  between  the  flow  and  combustion  is  also  represented  in  the  boundary  conditions  in  each 
flamelet  calculations.  The  coupling  in  the  other  direction  is  done  as  before,  i.e.  as  volumetric 
expansion  and  baroclinic  vorticity  sources  along  the  flame  front.  We  are  currently  implementing 
this  model  in  the  vortex  simulation  of  the  reacting  shear  layer  and  comparing  its  results  with  those 
obtained  using  the  transport  element  method  We  have  also  continued  our  investigation  on  the 
effects  of  density  variation,  and  its  coupling  with  forcing  on  the  growth  of  the  layer.  We  find  that 
while  the  growth  of  the  unforced  shear  layer  is  a  monotonic  function  of  the  density  ratio,  that  of  a 
forced  shear  layer  is  determined  by  the  momentum  ratio.  It  reaches  a  minimum  when  the  two 
streams  have  equal  momenta  and  grow  substantially  under  other  conditions. 
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I.  OBJECTIVE 

The  objectives  of  this  research  are: 

(1)  The  development  and  application  of  accurate  and  efficient  numerical  methods  for  the  integration 
of  the  time-dependent,  three-dimensional  Navier-Stokes  equations,  the  energy  and  species 
continuity  equations  at  high  Reynolds  and  Peclet  numbers,  naoderate  to  high  Damkohler 
numbers,  and  high  rates  of  heat  release.  Solution  methods  are  constructed  for  reacting  shear 
flows  at  moderate  to  high  Mach  numbers. 

(2)  The  investigation  of  mechanisms  of  flow-combustion  interactions  on  the  basis  of  the  results  of 
the  numerical  simulations,  and  the  study  of  how  these  interactions  can  be  exploited  to  control 
the  combustion  processes  in  turbulent  shear  flows.  Numerical  simulation  results  are  also  used 
to  construct  physical  models  which  can  be  used  in  turbulent  combustion  closure. 

Our  work  focuses  on  the  development  of  grid  free,  Lagrangian  field  methods.  In 
particular,  we  have  extended  the  application  of  the  vortex  element  method  and  formulated  the 
transport  element  method  for  the  simulation  of  compressible  reacting  flow.  We  have  also 
formulated  a  compatible  flame  sheet  model  for  the  simulation  of  high  Damkohler  number 
combustion.  To  validate  the  numerical  schemes  and  analyze  flow-combustion  interactions  in  shear 
flows,  we  have  focused  on  simulations  of  shear  layers  in  2D  and  3D,  with  particular  emphasis  on 
non-premixed,  variable  density  cases.  In  analyzing  flow-combustion  interactions,  we  obtain 
solutions  in  which  the  effects  of  the  velocity  gradient,  density  gradient,  heat  release  rate,  pressure 
gradient  and  confinement  are  varied  individually. 


II.  PERSONNEL 

During  the  funding  period  of  1990-1991,  two  graduate  students  and  a  postdoctor  were 
supponed  by  this  program: 

(1)  Marios  Soteriou  is  expected  to  finish  his  Ph.D.  thesis  by  the  end  of  this  academic  year  on  the 
Simulation  of  Combustion  in  a  Reacting  Shear  Layer. 

(2)  Van  Luu  started  his  Ph.D.  work  this  summer  on  the  Simulation  and  Mixing  and  Combustion 
in  a  Reacting  Shear  Layer  with  High  Rate  of  Heat  Release. 

(3)  Omar  Knio  spend  a  9-month  as  postdoctor  (50%  support  from  this  project)  extending  his 
work  on  the  three-dimensional  vonex  schemes  to  handle  compressible-shear  flows. 
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m.  WORK  STATUS 


ni.l.  Unsteady  Flame  Sheet  Model  for  Reacting  Flow 

The  simulation  of  combustion  in  the  reacting  shear  layer  has  so  far  been  accomplished 
using  the  transport  element  method  in  which  the  species  continuity  equations  are  integrated  along 
with  the  momentum  and  energy  equations  to  determine  the  local  rate  of  reaction  at  different  points 
in  the  field.  This  method  provides  a  direct  simulation  type  approach,  i.e.,  assumptions  regarding 
the  structure  of  the  reaction  zone  and  the  modes  of  interaction  between  the  flow  and  combustion  are 
not  made  apriori.  Instead,  the  simulations  have  been  used  to  determine  the  structure  of  this  zone 
under  conditions  of  low  and  moderate  Damkohler  numbers. 

As  the  Damkohler  number  increases,  it  has  been  found  that  the  thickness  of  the  reaction 
zone  becomes  much  smaller  than  the  typical  flow  scale  as  determined  by,  e.g.,  the  local  vorticity 
thickness  of  the  vorticity  layer.  This  substantially  increases  the  number  of  computational  elements 
required  to  resolve  both  the  flow  and  species  gradients  simultaneously  leading  to  a  sharp  rise  in 
computational  cost.  In  the  application  of  a  direct  method  such  as  the  transport  element  method  one 
needs  to  transpon  all  the  species.  This  increase  in  cost  is  compounded  by  the  need  to  extend  the 
models  to  accommodate  multi-step  chemical  kinetics.  If  a  flame  sheet  model  is  sufficient,  then 
only  the  local  concentrations  of  the  minor  species  are  needed. 

The  unsteady  flame  sheet  model  we  have  developed  will  be  described  in  detail  in  an 
upcoming  publication  (which  will  be  submitted  to  the  24th  Symposium).  It  is  based  on  assuming  a 
one-dimensional  flame  structure  in  the  direction  normal  to  the  local  flame  front  The  formulation  is 
developed  as  an  extension  of  the  Carrier,  Fendell  and  Marble  original  formulation  of  a  strained 
flame  with  infinite  rate  chemistry.  In  our  model  we  allow  for  finite  diffusion  and  finite  kinetic 
rates.  We  have  used  the  model  to  study  the  effect  of  unsteady,  primarily  periodic  strain  on  the 
flame  development  and  to  prove,  or  otherwise,  some  of  the  basic  assumptions  in  the  flame  sheet 
closure  models  of  turbulent  combustion.  The  results  are  summarized  next. 

In  the  flame  sheet  closure  model,  one  assumes  that  the  total  burning  rate  is  the  local  burning 
rate  ccnresponding  to  the  mean  strain  rate  multiplied  by  the  average  area  of  the  flame.  Using  the 
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mean  strain  to  compute  the  burning  rate  assumes  that  the  flame  does  not  respond  to  the  unsteady 
component  of  the  strain.  We  find  that  to  be  incorrect  in  cases  when  the  local  strain  exceeds  the 
quenching  strain  and  when  the  frequency  of  the  oscillating  flame  is  comparable  with  the  flow  time 
scale.  A  strain  generated  by  a  turbulent  flow  can  be  decomposed  into  two  components:  a  mean 
value,  Em,  and  a  root-mean-square,  e'.  When  the  flame  is  strained  such  that  the  maximum  strain, 
Em+e’  <  Eq,  where  Eq  is  the  mean  strain  which  quenches  a  flame,  we  find  that  applying  the  results 
of  steady  calculations  leads  to  a  reasonable  approximation  of  the  burning  rate.  However,  when 
Em+E'  >  Eq,  substantial  reduction  of  the  burning  rate  below  the  value  computed  using  the  mean 
strain  and  the  possibility  of  complete  extinction  are  encountered.  The  reduction  depends  on  the 
frequency  of  the  oscillations.  Moreover,  when  Em-E'  <  0,  the  reaction  zone  may  experience 
extinction  due  to  the  rapid  depletion  of  one  of  the  reactants.  These  results  are  very  important  in 
leaction  closure  models. 

We  are  currently  implementing  this  unsteady  flame  sheet  model  into  the  vortex 
computations  and  comparing  its  results  with  those  obtained  using  the  transport  element  method. 
Results  are  being  documented  in  a  paper  which  will  be  presented  at  the  Aerospace  Sciences 
Meeting  in  January,  1992. 

III.2.  Baroclinic  Effects  in  a  Variable-density.  Spatially  Developing  Laver 

We  continued  our  investigation  on  the  effects  of  variable  density  across  the  shear  layer  on 
its  growth  rate,  entrainment  ratio  and  motion  of  individual  eddies  within  the  mixing  zone.  During 
this  period,  we  focused  the  effect  of  upstream  forcing,  in  the  form  of  harmor.ic  modulation  of  the 
flow  rates,  <mi  these  parameters  by  comparing  results  of  an  unforced  layer  with  those  of  a  forced 
layer.  There  are  some  striking  qualitative  and  quantitative  differenres  between  the  two  cases.  The 
preliminary  results  have  been  included  in  recent  AIAA  paper,  with  is  attached  to  this  report.  The 
main  ccmclusions  are  listed  below. 

In  the  unforced  layer,  there  is  a  monotonic  dependence  of  the  growth  rate  on  the  density 
ratio  between  the  two  streams.  As  the  density  o^  the  slow  stream  increases  beyond  that  of  the  high 
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speed  stream,  the  growth  rate  of  the  layer  increases,  the  volumetric  entrainment  ratio  favors  the 
high  speed  fluid  and  individual  eddies  move,  with  respect  to  the  mean  flow  speed,.in  the  direction 
of  the  slow  stream.  These  results  are  in  agreement  with  the  experimental  measurements  and  agree 
also  with  the  conclusions  of  the  linear  stability  theory  of  variable-density-shear  layers.  Models 
based  on  the  redistribution  of  the  vorticity  within  the  large  structures  by  the  action  of  the  baroclinic 
torque  were  constructed  to  physically  explain  the  observed  trends,  which  as  mentioned  before  are 
well  correlated  by  the  density  ratio  (at  a  constant  velocity  ratio). 

In  the  forced  layer,  there  is  an  interesting  and  unexpected  non-monotonic  dependence  of  the 
growth  rate  on  the  velocity  ratio,  although  the  other  two  characteristic  parameters,  the  entrainment 
ratio  and  eddy  motion,  behave  in  the  same  way  as  in  the  unforced  case.  The  best  parameter  space 
which  we  found  for  characterizing  the  behavior  of  the  layer  in  this  case  is  the  momentum  ratio.  We 
found  that  the  smallest  growth  rate  (at  all  velocity  and  density  ratios  investigated)  occurs  at 
momentum  ratio  of  unity.  The  growth  rate  increases  for  momentum  ratios  smaller  and  larger  than 
unity.  This  behavior  has  been  observed  before  in  the  case  on  a  jet  issuing  into  an  atmosphere  with 
a  different  density.  The  mechanism  leading  to  this  critical  dependence  on  the  momentum  ratio  was 
found  to  be  the  pairing  of  the  eddies  which  can  no  Icwiger  be  described  as  vortex  monopole  but  as 
vortex  dipoles. 

Future  work.on  this  subject  will  consider  the  effect  of  different  forms  of  forcing  using 
different  combinations  of  the  fundamental  and  sut^armonic  frequencies  organized  at  different 
phasing  relationships.  Extention  to  3D  spatially  developing  layer  will  follow. 

ni.3.  Mixing  and  Combustion  between  two  Concentric  Jets 

With  partial  support  from  the  Gas  Research  Institute,  we  have  applied  the  methodologies 
developed  during  the  course  of  this  research  project  to  investigate  the  mechanism  of  mixing 
between  two  concentric  jets  and  the  formation  of  large  unsteady  structures  which  appear  to  govern 
the  burning  process.  This  problem  has  been  the  subject  of  several  experimental  investigations  in 
the  Air  Force  Laboratories  and  Sandia  Laboratories  since  it  represents  a  generic  configuration  of 
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combustors  in  propulsion  and  industry.  It  also  lends  itself  well  to  the  analysis  techniques  which 
we  have  developed. 

The  results  of  this  project  have  appeared  in  several  recent  publications.  One  of  these  papers 
is  attached  in  the  back  of  this  report.  We  find  that  as  the  velocity  of  the  inner  jet  increases,  it 
penetrates  the  recirculation  zone  of  the  outer  flow  and  causes  its  shear  layer  to  become  unstable. 
This  instability  is  manifested  in  the  shedding  of  large-scale  structures  form  the  shear  layer  between 
the  outer  recirculation  zone  and  the  free  stream.  As  the  jet  velocity  increases  further,  the  inner 
stream  breaks  through  the  recirculation  zone  and  large  eddies  are  shed  from  both  the  inner  and  the 
outer  shear  layers.  The  eddies  are  subsequently  paired,  forming  mixed  structures  downstream. 
These  structures,  which  form  at  a  finite  regular  frequency  that  depends  on  the  velocity  ratio  of  the 
two  streams,  resemble  the  flame  turbules  which  have  been  observed  in  the  reacting  flow 
experiments.  We  plan  to  extend  these  computations  by  adding  a  combustion  model  to  understand 
how  heat  release  affects  this  mixing  mechanism. 
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